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The missing metals problem. I. How many metals are in 
submm galaxies? 
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ABSTRACT 

At redshifts larger than 2, a large fraction (80%) of the metals are apparently yet 
undetected. We use a sample of sub-mm selected galaxies (SMGs) with molecular 
gas and dynamical mass measurements from the literature to put constraints on the 
contribution of such galaxies to the total metal budget. Compared to Lyman break 
galaxies (LBGs), for example, SMGs are rarer (by a factor of 10 or more), but contain 
much more gas and are more metal rich. For SMGs brighter than 3 mJy, we estimate 
that SMGs contain only <9% of the metals when we combine the observed dynamical 
masses (< Md yn >~fewxl0 n M Q ), number density (n ~ 10 -4 Mpc -3 ), observed gas 
metallicity (Z ~1-2Z & ), and observed gas fractions (/ gas ~ 40%) assuming a molecular 
to neutral hydrogen ratio of 1. Including SMGs fainter than 3 mJy, we estimate that 
SMGs contain about ^15% of the metals, where our incompleteness correction is 
estimated from the dust mass function. Our results are strong upper limits given that 
high gas fractions and high overall metallicity are mutually exclusive. In summary, 
SMGs make a significant contribution to the metal budget (<15%) but not sufficient 
to solve the 'missing metals problem'. A consequence of our results is that SMGs can 
only add rj 3.5% to JIdla, and can not be the source of a significant population of 
dusty DLAs. 
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1 INTRODUCTION 



A direct consequence of star-formation and, in particu- 
lar, of the star- formation history (SFH) (|Lillv^^alJ [T 996t 
iMadau et al.1 19961: iGiavalisco et alJ2004lHopkinsl2004 and 

others) is the production of heavy elements, known as met- 
als. Indeed, for a given initial mass function (IMF), the 
total expected amount of metals pz, e xpected formed by a 
given time t is simply the integral of the star formation 
density (SFD) (>*( €) times < p z >, where < p z > is th e 
mean stellar yield JSongaila et alJll990l:lMadau et alJll99rj) : 
pz.cxpoctod =< Pz > Jgdt p*(t) . Using a Salpeter IMF 
and the type II stellar yields ( for solar meta l licity) p z (m) 
from lWooslev fc Weaver! Il995t> . IMadau et alJ Jl996h found 
that < p z >— or 2.4%. Using th e SFH parameteri zed 
(in a LCDM cosmology) either as in ICole et ail (|200lT l or 
by a constant star formation rate (SFR) bey ond z = 2, we 
find that the total co-moving metal density is feouche et alJ 
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PZ, expected 



4.0 x 10 6 M Q Mpc" 3 , 



(1) 



after integrating the SFH over the redshift z range from 4 
to 2. This is about 25% of the 2 = metals. 

But at redshifts z ~ 3, our knowledge of the cosmic 
metal budget is still highly incomplete. Indeed, only a small 
fraction (20%) of the budget is actually seen when one adds 
the contribution of the Lya forest (A r jj I = 10 13-17 cm -2 ), 
damped Lya absorbers (DLAs) (Ajjj > 10 20 3 cm" 2 ), and 
galaxies such as Lyman break galaxies (LBGs) |pettini et alJ 
Il999t |PageJl2002l:lPettinil200a I Wolfe et alJl2003h . 

To account for the remaining 80% of Eq.0 or the "miss- 
ing metal^" there are two likely possibilities, as pointed out 
bv lPettinil fcOO&t . Either they are in a galaxy popula tion 
not yet accounted for in the budget of iPettini l2003Tl . or 
th ey are in a hot phase which is currently difficult to detect. 
In lBouche et alJ <l2005l . hereafter paper III), we discuss fur- 
ther the missing metal problem and th e latter alternative. In 
this paper and in iBouche et ail (120051) (paper II), we discuss 
the former. In paper II, we discuss the contribution of bot h 
the z ~ 2.2 UV selected galaxies, "BX" dSteidel et al.l2004h . 
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and near-IR selected galaxies (e.g. lFranx et al.l2003l) . In this 
paper, we discuss submm selected galaxies (SMGs). SMGs 
are potentially good candidates for hiding met als given that 
they a re both gas and metal rich. For instance. iDunne et alJ 
(2003, hereafter D03) explored the contribution of SMGs to 
the metal budget using the dust mass function (DMF) of 
high redshift (z ~ 2-3) submm galaxies constructed from 
deep blank field SCUBA surveys. From the DMF, D03 in- 
ferred, using chemical models, the co-moving density of met- 
als and baryons associated with the ISM of submm galaxies 
(SMGs). They concluded that all of the remaining metals 
(80%) are in the ISM of SMGs and that the mere existence 
of SMGs is enough to close the metal budget. 

Recently, the gas content and metallicity of SMGs 
have been directly est i mated or constrained in a few 
cases IGenzel et alJhoollNeri et all2003| : lGreve et alJ20ol 
iTecza et alJ 12004 ISwinbank et alJ 12004 . These measure- 
ments now allow us to put more direct limits on the contri- 
bution of SMGs to the cosmic metal budget. In this paper, 
we will show that indeed SMGs contribute significantly to 
the metal budget, but their contribution is <10% and, even 
optimistically, cannot be more than ~ 20%. In paper II, 
we show that z ~ 2.2 galaxies contribute significantly to 
the metal budget, up to 15-20%. Thus, combining all the 
know galaxy populations at z > 2.0, there is about 50% 
of the metal budget that has been accounted for. In pa- 
per III, we will explore whether the remaining metals have 
been expelled from small galaxies into the IGM (i.e. into a 
hot non-detectable phase). 

In section [3] we compare our results to those of D03 in 
a ACDM cosmology. 

In the remainder of this paper, we used Ho = 
70 h 70 kms" 1 Mpc" 1 , fi M = 0.3 and Q A = 0.7. 



2 HOW MANY METALS IN SMGS? 

In this section, we use the recent observa tions of gas content 
(H2) and dynamical ma sses of 7 SM Gs dGenzel et alJl2003t 
ITecza et alJl2004l: iGreve et alJl2005l) to put constraints on 
the contribution of the SMGs to the metal budget. We first 
summarize the observations. 



2.1 Properties of SMGs 

Recent measurement of the dynamical and gas masses 
of currently now a dozen z ~ 2.5 submm sources have 
been made from CO line emission (some resolved) using 
both OV RO and the IRAM Plateau de Bure i nterferom- 
eter (e.g. | Fraver et alJ|l998l : iGenzel et aT]|2003l : iNeri et alJ 
l2003l: lGreve et al.ll2005l) . Thev have redshifts spanning 1.0- 

1O 13 Ij0, and have large 



3.3, bolometric luminosities Lboi 

molecular mass M gas >2x 1O 1O M0 and dyna mical masses 
M d vn > 0.5 x 10 11 M Q . From the compilation bv lGreve et alJ 
(2005), we find that out of the dozen SMC with CO detec- 
tions, 7 meet the following two criteria: (i) 2 > 2 and (ii) an 
intrinsic (de-lensed) Ssso flux >3 mjy. The redshift cut-off 
is natural given the aims of this paper, and the flux thresh- 
old corresponds to the one used by D03 (see section |3J. We 
note that SMMJ14011+1152, with an intrinsic S 8 so flux of 
2.9 mjy, could be included in our sample if the magnifi- 



cation is slightly smaller. Our mean values do not change 
significantly if one includes SMM14011+1152. 

The properties of these sources are listed in Table 
One can see from the table that, on average, ~20% (and up 
to 50% for J01411) of the dynamical mass of SMGs is made 
of molecular (H2) gas. Excluding SMMJ01411, the averaged 
gas mass, velocity width and dynamical masses are M gaa = 
4.0xl0 10 M Q ,FWHM~ 700 kms" 1 , M dyn ~ 2.1 xlO 11 M Q , 
respectively (see the bottom of Table Q. 

Very few SMGs have had their gas phase metal- 
licities mea s ured or constrained. In the case of J14011, 
ITecza et alJ 120041) used the near-infrared integral field spec- 
trometer, SPIFFI (now SINFONI) on the ESO-VLT to 
measure the nebular emission line ratios of J 14011. Us- 
ing the classical optical diagnostic r atio, R23(=[OIIlAA 3726, 
3729 + [OIII]AA4959, 5007/H/3), ITecza et alJ <l2004h in- 
ferred a metallicity o f +0.27l°; 1 5dex (Z ~ 1.9±g;| Z Q ). 
ISwinbank et alJ (120041) used long-slit spectroscopy to mea- 
sure the [Nil] / Ha ratio to infer metallicities using the cali- 
bration of lPettini fc Paeel l2004h . The median of their sam- 
ple is slightly below solar. Broadly speaking, SMGs have 
metallicities close to solar and up to ~ 2 Zq . 

In order to estimate the total metal contribution of 
submm galaxies, it is necessary to estimate the true num- 
ber density of SMGs corrected for the "duty cycle," the 
fraction of cosmic time over which submm galaxies are ob- 
served. IGenzel et alJ (I2003F) estimated the raw co-moving 
density n of SMGs to be ~ 10" 5 h ro Mpc" 3 from the ob- 
served area covered by the SCLS and the estimated red- 
shift range 1 < z < 5. The duty cycle of SMGs can 
be constrained dire ctly from the gas depletion time scale. 
IGenzel et al.l (J^003) showed that such a luminous galaxy 
(with a SFR ~ SOOMoyr" 1 ) would use up its gas in ap- 
proximately 4 x 10 8 yr and make ~ 2 x 10 11 Mq of stars. 
ITecza et all <2004h estimated an age for J1401 1+0252 of 
^200 Myrs from the strength of its Balmer break. Us- 
ing a SMG time scale of 4 x 10 8 yr, and assuming that 
SMGs are evenly distributed over l^z^5 (~ 4.5 Gyr), this 
would imply that SMGs are 'on' 10% of the time or have a 
duty cycle r of 0.1. Thus, the co-moving density is at least 
n > 1.3 X 10~ 4 fe?n f^)" 1 M P c " 3 t iGenzel et al]l2003h . 
Independently, IChapman et al.l 120051) estimated a raw co- 
moving density n of z ~ 2.5 SMGs above Lfir = 10 12 ' 5 
of n ~ 10~ 5 /i 3 o Mpc~ 3 using m e asured redshifts of radio 
identified SMGs. IChapman et ail 120051) modeled the red- 
shift distribution of SMGs with a Gaussian distribution of 
width a z ~ 1.3 (covering ~ 1 Gyr), and assumed a time 
scale of 10 8 yr, yielding a similar duty cycle r of 0.1. Thus, 
the true co-moving space density estimated by two groups, 
n=l-3xHT 4 h 70 (f^-)" 1 Mpc" 3 . 

Another way to estimate the duty cycle of SMGs 
is to use their clust ering strength (as pointed out by 
IChapman et alJ l2005l0 and Fig. Q Both panels show 
the comoving number density of sources vs. redshift. In 
the left panel, lines of constant halo mass log Mh = 
14, 13, 11, 10, 9 are shown. In the right panel, lines of 
constant clustering amplitude ro = 15, 10, 8, 6 , 4, 2 Mpc 
(from bottom to top) are shown. iBlain et al.l i2004h es- 
timated the SMG auto-correlation length ro to be ro = 
7 h^Q Mpc (represented by the star), from which one would 
infer a halo co-moving abundance density of n ~ 1 x 
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Table 1. Gas masses and dynamical masses o f S MGs with z > 2 an d a n intrinsic flux 5850 gre ater than 3 mjy. T he a verage quantities 
are sh own. References: (Fl lOenzel et alJ 120031) . f2^ lNeri et aflfcOOaft . CS lGreve et alJ bOO-fl . f4l lTecza et alJj2004) . (5l lSwinbank et alJ 
l2004\ . 



Name S 850 a zco M gas b FWHM c M dyn d Refs [O/H] e Refs 

(mjy) x!O lo M kms" 1 xlO lo M 



SMMJ02399-0136 


9.6 


2.8076 


6.0 


1100 


60 


1 






SMMJ09431+4700 


8.8 


3.3460 


2.0 


420 


5/7 


2,3 






SMMJ131201+4242 


6.2 


3.408 


4.2 


530 


12 


3 






(SMMJ14011+1152 * 


2.9 


2.5652 


3.4 


190 


6 


1,3 


0.3 


4) 


SMMJ16358+4057 


8.2 


2.3853 


5.6 


840 


9/35 


2,3 






SMMJ16366+4105 


10.7 


2.450 


4.6 


870 


9/37 


2,3 


0.1 


5 


SMMJ16371+4053 


10.5 


2.380 


2.4 


830 


34 


3 


-0.1 


5 


SMMJ22174+0015 


6.3 


3.099 


3 


780 


28 


3 






average f 


8.6 




3.97 ± 1.55 


767±226 


21±21/30±17 









"Intrinsic submm fluxes. 
^Molecular gas masses from CO line emission. 
c Full width at half maximum of the CO line. 
^Dynamical masses for the inclination i = 45 and using /170 = 1. 
e [0/H] metallicities using (O/H) from lAsplund et alJ i2004h . 
■^The amplification of this source was revised from 2.5 to ~5 JSmail fc et aljfeoOSl) moving it just below our threshold of 3 mjy. 

^Excluding SMM14011+1152. 



fi =0.3, A=0.7, h = 0.7, ct =0.8 




z z 



Figure 1. Left: Co-moving number density of sources vs. redshift. Lines of constant halo mass are shown by the contours for logMj, = 
14, 13, 11, 10, 9, from bottom to top. The filled square shows the observed number density of SMGs, n ~ 1 X 10" 5 h 3 Mpc — 3 , with the 
uncertainty and redshift coverage represented by the dark shaded area. Right: Same as left, but the contours represen t lines of constant 
clustering amplitude ro = 15, 10, 8, 6, 4, 2 Mpc (bottom to top). The correlation length of SMGs (ro ~ 7 Mpc lBlain et ahlEoolF) 
is shown by the star. The halo mass inferred from the clustering is ~ 10 12 5 Mq. Given that the observed n (left panel, solid square) is 
lower than the one inferred from the clustering (star) by a factor of ~ 10, the duty cycle of SMGs is about ~ 0.1. The number density of 
SMGs correc ted for the du t y cycl e is shown by the light shaded area in both panels. Both panels were produced using the Press-Schechter 
formalism of lMo fc White! fcoOSft . 



1CF 4 /170 Mpc" 3 (light shaded area). This is «10 times larger 
than the observed number density (n ~ 1 x 1CP 5 Mpc" 3 , 
dark shaded area). These two numbers can be reconciled if 
SMGs are short lived, with a "duty cycle" of the order of 
0.1. Nonetheless, the several methods of estimating the duty 
cycle of SMGs agree to within a factor of a few. 



2.2 Consequences of SMGs properties 

From the observed properties of SMGs (summarized in Ta- 
ble0, in this section, we derive the co-moving baryonic den- 
sity and co-moving metal densities in SMGs (summarized in 
Table 0. 

The dynamical masses of SMGs cover the range 2 — 
3 x 10 11 ftj Mg (Table and assuming that this mass is 
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entirely baryonic, the co-moving baryonic density in submm 
galaxies is observed to be: 

p b , ohs ~ 2.7 x 10 7 ^ot^-r 1 M Q Mpc" 3 . (2) 

The observed mean gas mass M gas of ~ 4.0 x 10 10 M 
corresponds to a gas fraction f g = 20% and implies a 
gas (molecular) co-moving density of p gas ,SMG = 5.2 x 
10 6 /i? o M Mpc -3 . 

Naturally, the ISM of the SMGs contains also an un- 
known amount of neutral gas (Hi). Since we do not know 
the neutral to molecular ratios, it is difficult to make a robust 
estimate of the total potential res e rvoir of gaseous material. 
In the local Universe. iKeres et alJ (2003) have estimated the 
co- moving mass contribution of the cold molecular gas (H2) 
and compared it to the atomic phase (H 1) of the interstellar 
media of galaxies. They find an average ratio of H2 to Hi of 
about 0.5 to 0.65. The molecular-to-neutral ratio increases 
with increasing H2 mass and is about 2 for galaxies with 

M H 2 ~ 109 M °- At 

very high infrared luminosities, similar 
to what is found for the SMGs, galaxies show H 2 to H I mass 
ratios of ~ 4 to 20 llMirabel fc Sanderdll989h . Using these 
local results, in order to provide a more robust estimate of 
the total metal content, we assume that (i) there is as much 
HI as molecular gas, which is likely to be an u pper limit 
given the local ratios iMirabel fc Sandersl lT989). and that 
(ii) the metallicity of SMGs is on average < Z >~ 1.9 Z@, 
the metallicity of SMM14011+1152. We find that the metal 
co-moving density in SMGs is: 

pz,obs < 3.6xl0 5 h 2 70 ^ M Mpc' 3 ,(3) 

or < 9% of the metal budget. 

We also note that our assumption that there is as much 
HI as molecular gas implies that SMGs can only add « 3.5% 
to the neutral gas content of the univer se Odla, and thus 
SMGs can not harbor the dusty DLAs of lVladilo fc Perouxl 
(2004). The two are distinct populations. 

It is absolutely necessary to set a firm upper limit to the 
contribution of SMGs to the cosmic metal density. This is es- 
pecially true in light of the claim by D03 that the metal con- 
tent of SMGs is sufficient to solve the missing metals prob- 
lem (more on this in §EJ. Eq. |3]is a hard upper limit from 
two independent lines of arguments. First, it is the maxi- 
mum gas-mass fraction allowed by the metallicity. Indeed, 
high metallicity and large gas frac tions are m utually exclu- 
sive in chemical evolution models. Ed mund's! (ll99CT) showed 
that when the gas fraction is much larger than 50% solar 
metallicity canno t be reached. In add i tion, the chemical evo- 
lution models of lEdmunds fc Ealesl <1998ft (used by D03), 
showed that the stellar contribution to the metal mass bud- 
get of SMG is small (approximately 1 /4 of Eq. |3J given the 
lower metallicity of the stellar component compared to the 
gas component (see section 3). Second, our mean metallic- 
ity of < Z >= 1.9 is large. If we assume Z = Zq and treat 
both the gas and the stellar component equally, the SMG 
contribution to the metal budget would remain the same. In 
this case, Eq. |3] is a strong upper limit given that 100% of 
the dynamical mass Md yn is in baryons (i.e. neglecting any 
contribution from dark matter to the dynamical mass). 

Thus, we conclude that Eq. |3 is a strong upper limit 
given our assumptions. SMGs with Ssso > 3mJy can not 
contribute more than ~9% of the z = 2 metal budget. 



3 THE DUST MASS FUNCTION AND ITS 
IMPLICATIONS 

D03 took a very different approach. They used submm data 
from published deep blank field SCUBA surveys and con- 
structed the dust mass function (DMF) of high redshift 
galaxies. From the DMF, they inferred, using chemical mod- 
els, the co-moving density of metals and baryons associated 
with the ISM of gala xies. They assum ed a dust tempera- 
ture of ~ 25 K from iPei et all |l999) . Because the dust 
mass is a strong function of the dust temperature, a lower 
(higher) dust temperature 20 K (30 K) will increase (de- 
crease) their mass estimates by a factor of about two (D03). 
D03 argued that most of the dust mass in SMGs would be a 
low temperature (Td ~ 20 K) component. We note that the 
recent observations of IChapman et all (120051) favor higher 
dust temperature (Td ~ 35 K). D03 also used the dust mass 
opacity measured at 125pm and extrapolated it to submm 
wavelengths using a depende nce with j3 = 2 (from their 
local survey iDunne fc EalesboOll) . We now discuss their re- 
sults, which are summarized in Table [2] 

3.1 DMF all 

D03 find that the DMF function is well described by a 
Schechter function with the parameters — 4.7 x 1O 8 M0, 
4>* d — 8.9 x 10 -4 Mpc -3 , and 01 = —1.08 in a A cosmology. 
The integral of the DMF gives the co-moving density of dust 
p d = T(2 + a)M2<j>d = 4.39 x 10 5 M Q Mpc" 3 x . From p d , 
they assumed that 40% (rj — 0.4) of the ISM metals are 
locked into dust grains, yielding a metal density: 

Pz.dmf = 1.1 x 10 6 h 70 — M Mpc" 3 , (4) 
V 

i.e., about 27% of the cosmic metal budget. The apparent 
discrepancy between that number and the original conclu- 
sion of D03, namely that SMGs contain > 70% of the cosmic 
metal density, is due to the different cosmology assumed. 

A strong lower limit comes from n — 1, i.e. assuming 
100% of the metals are locked onto the submm emitting dust 
grains. In that case, at least 11% of the metals are in the 
ISM of SMGs. Note that (i) this does not include the metals 
in stars, and (ii) it will not depend on the chemical evolution 
models discussed below. 

D03 used a "closed -box" chemical evolution model 
fedmunds fc Ealesl Il99gl) to convert the dust content into 
a total baryonic content. They find that SMGs contribute a 
baryonic co-moving density of about: 

p 6 ,DMF ^ 7.2 x 10 7 hro M Mpc" 3 . (5) 

In these chemical evolution models, they have assumed 
that the submm sources are observed at their maximum dust 
mass, i.e., at the peak of the dust mass to baryonic mass 
ratio. Thus their baryonic mass density is a strong lower 
limit. This peak occurs at a gas fraction of roughly 50%. 

1 D03 used a Qjvf = 1 cosmology throughout their paper. In the 
remainder of this paper, we used their DMF in a Qm = 0.3, 
A = 0.7 cosmology. This lowers their cosmic densities (pi, dmf 

.;S x 10 s ft.75 M Mpc" 3 , pz,dmf = 1.9 x 10 6 h 75 M Q Mpc~ 3 ) 
by a factor 1.7. The global factor between our cosmology and D03 
is 1.8 including the change from /17s to h 7 Q. 
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Table 2. Baryons and metal cosmic densities in SMGs from this paper. 



Baryons 


p(M Mpc" 3 ) 
5.98 x 10 9 h\ 


P/Pc 

0.044 ft" 2 


p/Pb (%) 
100 


Note 


Z/Zq 


Pz(M Mpc" 3 ) 


Pz/Pc 


Pz/pZ,tot (%) 


Stars 2 < z < 4 


1.21 x 10 s ft£> 


0.000890 h n f 


2.02 


/ SFR a 


1.25 b 


4.00 x 10 6 


2.94 x 10" 5 


100 


SMGs > 3mJy: 


SMGs Baryons 
SMGs ISM (H2) 
> ISM (H2+HI) 


2.73 x 10 7 h% 
5.16 x 10 6 ft 2 
1.03 x 10 7 h\ 


0.000201 h% 
0.000038 ft° 
0.000076 h% 


0.46 
0.09 
0.17 


fa = !9% 
fa = 38% 


1.9 


3.63 x 10 s 


2.67 x 10" 6 


9.1 



"Cosmic stellar density calculated from the integrated SFH taking into account a recycled fraction of R = 0.28. f g is the gas fraction. 
6 Averaged yield (= 1/42 = 1.25 Z ) for type II SN with m > 10 M iMadau et alll996ft . 



Table 3. Results from the dust mass function (D03). Numbers in bold are taken from D03 and corrected for our cosmology. 



Baryons 


p(M Mpc -3 ) 
5.98 x 10 9 h\ 


P/Pc 

0.044 ft" 2 


p/pb (%) Note 
100 


Z/Z 


p z (M Mpc~ 3 ) 


Pz/Pc 




Pz/PZ, tot (%) 


Stars 2 < z < 4 


1.21 X 10 s h^ Q 


0.000890 ft" 2 


2.02 J SFR a 


1.25 b 


4.00 x 10 6 


2.94 x 10" 


-5 


100 


SMGs dust (DMF) 
SMG baryons 

- ISM 

- Stars 


4.39 x 10 5 h£ 
7.65 x 10 7 h£ 

3.82 x 10 7 h\ 
3.82 x 10 7 h\ Q 


0.000003 
0.000562 ft" 1 
0.000281 h n l 
0.000281 


0.01 
1.28 

0.64 T) = 0.4 
0.64 


1.52 
0.33 


1.10 x 10 6 
2.41 x 10 s 


8.07 x 10- 
1.78 x 10" 


-6 
-6 


27.4 
6.0 


SMG bar. ( > 3mjy): 

- ISM 

— Stars 


4.47 x 10 7 ft 70 
2.24 X 10 7 ft 70 
2.24 X 10 7 h\ 


0.000329 ft" 1 
0.000164 h~l 
0.000164 ft" 1 


0.75 
0.37 
0.37 


- >0.33 


6.42 x 10 5 
1.41 x 10 5 


4.72 x 10" 
1.04 x 10" 


-6 
-6 


16.0 
3.5 



"Cosmic stellar density calculated from the integrated SFH taking into account a recycled fraction of R = 0.28. ij = 0.4 is the fraction of 

the ISM metals that are assumed to be locked onto dust grains. 
b Averaged yield (= 1/42 = 1.25 Z e ) for type II SN with m > 10 M iMadau et alll996ft . 



As a side note, the inferred metallicity from D03 is close 
to the observed metallicity of SMGs. From the metallicity 
density (Eq. |1J and the baryonic density (Eq. |KJ obtained 
from the chemical models) , the mean metallicity of the ISM 
is pz,DMF/{Pb,DMF • fg) or about ~ 1.51 Zq, assuming that 
all the metals are in the ISM and a 50% gas fraction f g . This 
mean metallicity is very close to the observed metallicity of 
one SMG (J14011) discussed in section 

3.2 DMF bright 

In order to compare this prediction to the observed prop- 
erties of SMGs discussed in section [5] one needs to com- 
pare the cosmic baryonic and metal densities of flux se- 
lected SMGs to the DMF with a similar flux limit. D03 
integrated the DMF down to Ssso > 3 mjy, but quoted only 
the stellar mass density (p s t ar = 2.24 x 10 7 ft.70 M Mpc -3 ; 
half of the baryons) and stellar metal density (pz,stm = 
1.4 x 10 5 h 70 M Mpc" 3 ) 2 . Thus, the cosmic baryon den- 
sity for bright (with Sgso > 3 mjy) submm sources of (twice 
the stellar density): 

2 We again converted their number (p;, dmf = 2 X p* = 2 x 3.8 X 
10 7 h 70 M Mpc" 3 ) in a CIm = 1 to a A cosmology. 



Pb,DMF,3mJy ^ 4.2 X Iff ft 7 oM M P C" 3 (6) 

i.e., ~ 0.6 times the number quoted in Eq. providing the 
completeness factor. We then scale the total metal cosmic 
density (Eq. [1| by this 0.6 factor to infer the metal cosmic 
density in sources with Ssso > 3m Jy: 

PZ,DMF,3mJy ~ 6.4 x 10 s ft 7 oM Mpc -3 , (7) 

or about ~16% of the estimated cosmic metal density 
(Eq.[TJ. 

If one compares the baryonic and metal densities pre- 
dicted from the DMF (Eqs.|S||7|l with the observed baryonic 
and metal densities (Eqs.|5J|3J, we conclude that the predic- 
tions from the DMF (in a A cosmology) were over-estimated 
by a factor of at least 2. This factor can be easily accounted 
for if one uses a higher dust te mperature (Td — 35 K ), as 
indicated by the observations of IChapman et al] (I2005T) . 



4 SUMMARY & DISCUSSION 

SMGs are gas- (gas fraction 20-50%) and metal-rich galaxies 
(Z/Zq > 1). Therefore, they are potentially good candidates 
for harboring the missing metals. From the observed gas 
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fraction ~ 40%, dyna mical mass Md vn a nd mean metallicity 
of > Z Q (supported bv lTecza et all2004l) of 7 SMGs brighter 
than £850 > 3 mjy, we show that 

• Based on the dynamical masses of SMGs an assuming 
a 100% baryon fraction with < Z >= Zq, SMGs can not 
contribute more than 9% of the expected cosmic metal den- 
sity; 

• Based the observed high gas fractions and observed high 
ISM metallicities Z » Zq, SMGs can not contribute more 
than 9% of the expected cosmic metal density; 

• the total contribution of SMGs, correcting for incom- 
pleteness (section 3.2), is times the contribution of SMGs 
brighter than 3 mjy, or < 15%; 

• our results imply that SMGs can only add « 3.5% 
to OpLA- Thus SMGs cannot harbor the dusty DLAs of 
IVladilo fc Perouxl 12004^1 . 

Early estimates of the contribution of the SMGs to the 
metal budget from the DMF were overestimated. The dis- 
crepancy however is mainly due to the assumed cosmology 
and to the low dust temperature used by D03 (Td ~ 25 K). 
We do agree with the conclusion of D03 that SMGs make 
a significant contribution to the cosmic metal budget, just 
not enough to solve the "missing metals problem" . 

We are still far from closing the metal budget, however. 
In addition to <9% of the m etals that are in SMGs, 5% are 
in z ~ 3 LBGs dPettinill2003h . ~ 15-20% inz~ 2.2 galaxies 
(see paper II), 8% in the forest (but see paper III), and 5% 
m DLAs 3 . Taking our results on SMGs at face-value, and 
ignoring the issue of double-counting, roughly 50% of the 
metals have been accounted for (see also paper III). 

We are exploring two main ave nues in t rying to close the 
missing metals problem. Following iPettinll i2003T) . either an- 
other population of galaxies has not yet been accounted for 
or there is a significant reservoir of metals in the IGM that 
has not been detected. A substantial fraction of the miss- 
ing metals may be hidden in a very hot, collisionally ionized 
gas. Based on simple order-of-magnitude calculations, in pa- 
per III, we will discuss the possibility that the remaining 
missing metals could have been ejected from small galax- 
ies via galactic outflows into the IGM in a hot (T> 10 6 K) 
that is difficult to detect using observed properties of local 
galaxies. 



ACKNOWLEDGMENTS 

We thank T. Greve for providing the Ssso fluxes of the 
SMGs; A. Baker, L. Tacconi, D. Lutz and M. Pettini for 
helpful comments and discussions. We thank the referee for 
his/her comments that improved the quality of the paper. 



3 Dusty DLAs iVladilo fc Perouxll2004l) . missed in current spec- 
troscopic DLA surveys (cur rent DLA samples show small molec- 
ular and dust contents, e.g. lEllison et afll200lt iMurphv fc Liskd 
|2004|) could amount to an additional 17% (paper III). These would 
be in a separate population from SMGs given that the amount of 
Hi in SMGs is less than 3-4% of C D LA- 
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